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Three-dimensional non-Abehan anyons have been theoretically proposed to exist in heterostruc- 
tures composed of type II superconductors and topological insulators. We use realistic material 
parameters for a device derived from Bi2Se3 to quantitatively predict the temperature and mag- 
netic field regimes where an experiment might detect the presence of these exotic states by means 
of a cooling effect. Within the appropriate parameter regime, an adiabatic increase of the magnetic 
field will result in a decrease of system temperature when anyons are present. If anyons are not 
present, the same experiment will result in heating. 
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I. INTRODUCTION 

Majorana fermions (which arc responsible for the non- 
Abelian content of Ising anyonsi) have been proposed 
to exist in a number of condensed matter settings. One 
particularly intriguing development involves the Majo- 
rana states trapped within the hedgehog defects that 
emerge at the interface between a topological insulator 
and the mouths of vortices in a type II superconductor—. 
Since the realization of anyon statistics is generally be- 
lieved to be impossible beyond two dimensions, the nov- 
elty of this setup is its three dimensional (3D) nature. 
The paradox resolves upon realizing that the anyons at 
hedgehogs are not merely point-like objects, but singu- 
larities of a U{n)/0{n) order parameter which may be 
modeled by ribbons joining hedgehogs^. The ribbons 
contain additional twist degrees of freedom relevant to 
particle exchange: twisting a ribbon transforms the state. 
The anyons in the Teo-Kane model provide a projective 
representation of the ribbon permutation group^. Inter- 
estingly, this remnant of braiding in higher dimensions 
seems to be special to Ising anyons: Fibonacci anyons, 
for example, seem to have no 3D version. While these 
mathematical details are interesting, they arc subordi- 
nate to the question relevant to physicists: Do 3D anyons 
really exist or are they merely theoretical constructs? To 
the best of our knowledge, this paper provides the first 
proposal to settle this question in the lab rather than via 
theoretical arguments. 

Majorana fermions have proven challenging to unam- 
biguously detect in any dimension not to mention the 
exotic 3D context we have in mind here. For exam- 
ple, within the setting of quantum Hall systems, relevant 
to 2D non-Abelian anyons, experiments usually involve 
edge state transporl^i. Recently, however, several bulk 
probes have been proposed^i^. In particular, one sugges- 
tion exploits the huge ground state degeneracy inherent 
to the non-Abelian anyon system to produce a cooling 
effecli^. This can be understood in analogy to adiabatic 
cooling via spin demagnetization, but here the entropy 
reservoir is an anyon system rather than a spin system. 



Crucially, the non-Abelian entropy is temperature inde- 
pendent but proportional to the number of non-Abelian 
anyons, which, in turn, can be adiabatically (or, more 
precisely, isentropically) controlled by a magnetic field. 
Thus, an adiabatic increase in the non-Abelian entropy 
necessarily implies a decrease in the rest of the system's 
entropy. This is accomplished by a lowering of the sys- 
tem's temperature. While this idea has been explored 
in a 2D setting^, we expect it to work even better in 
3D where the topological and conventional sources of en- 
tropy (which are fundamentally 3D) will be on an equal 
footing. 

In this paper we explore the possibility of detecting 3D 
anyons via cooling by applying the Teo-Kane model to a 
specific heterostructure involving a topological insulator 
and a superconductor. Since we arc mainly interested in 
demonstrating feasibility, our principal goal will be to es- 
tablish an experimentally accessible temperature window 
Tl < T < Tjj in which the cooling effect is significant. 
We make quantitative estimates using realistic material 
parameters. The upper bound of the temperature win- 
dow, Tu, is determined by the temperature where sources 
of entropy other than anyons begin to dominate the sys- 
tem. The lower bound of the temperature window, Tl, 
is established by calculating the ground state energy de- 
generacy splitting that results from Majorana-Majorana 
tunneling. Operating at temperatures above this split- 
ting energy allows us to treat the low energy Hilbert space 
as essentially degenerate. 

In Section |TT] we briefiy mention existing proposals to 
detect the more conventional type of non-abelian Ising 
anyons below three dimensions. Section [IIII describes the 
physical structure and material components of the device 
for our 3D anyon system. The model we use to analyze 
this device is explained in Section ITVl from which we con- 
struct Majorana solutions in Section IVl and calculate the 
ground state degeneracy splitting as a function of Majo- 
rana separation in Section IVII Section I VIII describes the 
calculation of all contributions to the system entropy. We 
use this in Section IVIIII to determine the feasibility of de- 
tection by cooling and establish the parameter regime in 
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which to carry out the search for 3D non-abehan anyons. 
Our conclusions are summarized in IIXI 



II. MAJORANA DETECTION PROPOSALS 

FOR D < 3 

There are now many proposals for condensed mat- 
ter realizations of Majoranas in two dimensional sys- 
tems. Among them, the most actively studied thus far 
is the quantum Hall state at filling factor 5/2 where de- 
tection schemes have been proposed using both edgei^ 
and bulk^i^ probes. Some tantalizing evidence has 
been reported based on the former— liS. Other systems 
in which Majoranas may exist and detection strate- 
gies have been theoretically proposed include strontium 
ruthenatesii, helium-3i^, topological insulator based 
heterostructureai^, semiconductor heterostructures^^^, 
cold atoms systemai^, and quantum wire nctworksii"— . 
The experimental confirmation of the existence of these 
2D Majoranas has lagged somewhat the large number of 
theoretical propositions. 

Nonetheless, all these Majorana fermions are of the two 
dimensional SU(2)2 Ising anyon variety^. They are repre- 
sentations of the braid group. The possibility for anyons 
beyond two dimensions came as a bit of a shock^i^ 
due to certain long-standing trusted mathematical argu- 
ments. It turns out that no mathematical theorems need 
to be corrected because these 3D anyons do not provide 
representations of the braid group but rather provide pro- 
jective representations of the ribbon permutation group*^. 

We also mention in passing that the quantum wire net- 
works involve ID Majorana fermions in a sense, but their 
braiding (for example using T-junctions as proposed by 
Alicea et a.&-) still requires the spanning of 2D real space 
by the network. Furthermore, these ID Majoranas are 
still linked to the 2D anyon statistics tied up with the 
braid group rather than the 3D projective ribbon per- 
mutation statistics of 3D anyons. 

III. DEVICE STRUCTURE AND MATERIALS 

Consider an s-wave, type II superconductor sand- 
wiched between layers of a 3D topological insulator. The 
superconductor occupies the region — < z < with 
flanking layers of topological insulator for ^ < \z\ < L. 
This structure is repeated along the z-direction to build 
a superlattice and thus a true bulk 3D system. See fig- 
ured Due to this periodicity, our calculations will only 
need to consider a single superlayer. For the topologi- 
cal insulator, we choose Bi2Se3 appropriately p-dopcd so 
that the Fermi level (ep) is pushed down into the bulk 
gap. For the superconductor, we choose the n-doped ma- 
terial Cuo.i2Bi2Se3 which is known to be a strongly type 
II (ac = A/^ « 50) bulk superconductor—. The precise 
nature of the superconductivity in this material is not 
yet known, but a theoretical proposal has suggested that 




FIG. 1: (color online) Two superlayers of the heterostruc- 
ture consisting of alternating layers of topological insulator 
(TI) and superconductor (SC). The external magnetic field 
induces vortices (green tubes) in the SC layers that terminate 
at the interface with the TI layers resulting in 3D anyons 
(orange balls). In addition to the obvious three dimensional 
distribution of anyons in real space, these anyons provide a 
projective representation of the ribbon permutation group^ 
which further distinguishes them from the more well-known 
two dimensional anyon system (which forms a representation 
of the braid group). 



it may be unconventional^i. In the absence of further 
experimental data, we will assume it to be s-wave; this 
and other caveats are discussed further at the end of the 
paper. 

Existing experimental results on Cuo.i2Bi2Se3 (almost) 
supply the required minimal input needed for our calcu- 
lation. ARPES provides numbers for the Fermi energy, 
velocity, wavevector, and the spin-orbit gap2^: ep = 
400 mcV, hvF = 3.8 eV A, kp = 0.1 A~^ w ep/fivp, 
and Ago = 150 mcV. We also need the superconduct- 
ing gap (A) and coherence length (^). Based on the 
existing experimental data, there are two ways of de- 
termining A and ^. Both quantities could be inferred 
from the known superconducting transition temperature, 
Tc = 3.8 K23^ using the BCS relations A = 1.764fcBTc 
and ^ = Hvp/ttA. On the other hand, we could use 
the measurement of the upper critical field, Bc2 = 1.7 
Tesla— , to give a more direct estimate of the coherence 

length ^ = ^ 27^5*2 ^ ^'^'^ ^ independent of BCS the- 
ory. This could then be combined with the BCS relation 
A = hvp/n^ to give a value for the superconducting gap. 
The fact that these two methods do not agreo^i^ pro- 
vides a puzzle for the community. What is needed is a 
direct measurement of A. which has not been reported 
yet. Since this paper is chiefiy concerned with vortices we 
will trust the i3c2-derived value of ^, which is 139 A. To 
be consistent, this value of ^ is then used to set A = 8.7 
meV. 

Very recently, another indirect estimate of the gap has 
been made using specific heat data^^. Weak-coupling 
BCS theory does not fit the data, while a fit to strong- 
coupling BCS theory produces parameter values incon- 
sistent with other measures. Therefore, the main con- 
clusions to be drawn from this experiment are the bulk 
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nature of the superconductivity, the lack of nodes in the 
gap, and the fact that the gap is larger than would be ex- 
pected from BCS theory using Tc- We therefore continue 
to use Bc2 to estimate the value of the superconducting 
gap. 



IV. THEORETICAL MODEL 

The low-energy theory is an eight-band Dirac 
mode^^Sl^: H = \ j' (fx^'f'H'i' where 



n 



Hd a 

A* --Hd 



(1) 



with diagonal terms given by Hjj = cx ■ p — ep - 
i'^^/3Aso- We use the standard Dirac-Pauli representa- 
(t\ f t \ . f 1 



tion: a = 



The fermion is an eight-component object deriving from 
spin, orbital, and particle-hole degrees of freedom: = 
(-0''^, — iV'"^a2) = {'u\,U2,vl,vl), where each and Vi 
has two components. Particle-hole symmetry enforces 
Vi = ia -iu* . We have written the model in the same way 
as Ref. 1271 . which is unitarily equivalent to the models 
used in Refs.lH and[28l Variants of this model have been 
studied for many years^^. 

We use a Dirac Bogoliubov-de-Gennes (BdG) model 
for this material. The construction of the effective model 
for Bi2Se3 has been discussed by several groups, see for 
example^ and references therein. While higher momen- 
tum terms can sometimes lead to interesting physics, the 
situation we are considering only requires use of a Dirac 
model. 



MAJORANA SOLUTIONS 



and solve for the zero energy solutions of the BdG equa- 
tion H^' = 0. At the kink (z = -L/2) this leads to the 
following Majorana zero-mode wavefunctions^LS^ 







V J(«+i)/2(fcFp)e-/4e'("+i)''/2 ; 



(2) 



while at the upper interface {z = L/2) we have an anti- 
kink given by cyclically permuting the right hand side 
two steps. The expressions for Vi follow from particle- 
hole symmetry and A/" is a normalization constant given 



by = 



where K and E 



are the complete elliptic integrals of the first and second 
kinds. These expressions are valid inside the supercon- 
ductor. We ignore the exponential tail decaying into the 
topological insulator which is miniscule due to the large 
insulating gap. 

An exact treatment of this problem would numerically 
solve for the order parameter profiles self-consistently. 
See, for example, Ref. [U However, since we de- 
sire analytical expressions we make the following stan- 

j J • TC ^- -IT— f '^(p')dp' -jr^ f A,a(z')dz' 

dard simpufication: e ' e ""f - ^ ' = 



sech(/9/7r^)scch(z/^2) where p = y/x^ + y^ and = 
hvp/Aso ~ 25.3 A. Thus, two experimentally deter- 
mined parameters infiuence the localization of the Ma- 
jorana wavefunction: ^ and A^o (or, equivalently, as 
defined above). determines the localization in the z- 
direction while ^ sets the decay length in the xy-plane. 
Since ^ ^ ^z, we might already speculate that the in- 
plane coupling between Majoranas will be much more 
important than the coupling in the z-direction for the 
degeneracy splitting. We turn to this issue next. 



We combine the superconducting and spin-orbit gaps 
into a single 3-componcnt order parameter—: n = 
(RcA, ImA, Aso). Within the superconductor, a vortex 
excitation along the z-axis is achieved by imposing the 
following profile on the superconducting order param- 
eter: A{p,0,z) = A(p)e™^e(L/2 - \z\). The interface 
between the topological insulator and superconductor is 
specified by imposing the following kink profile on the 
spin-orbit gap: Aso{p,0,z) = +Aso for \z\ <^ L/2 and 
Aso{p, d, z) ~ —Aso for \z\ ^ L/2. In this way the band 
is inverted within the topological insulator, while taking 
the opposite sign in the (topologically trivial) supercon- 
ductor. Note that we have considered a topological insu- 
lator and trivial superconductor, but the kink would also 
exist at the interface of a trivial insulator and a topolog- 
ical superconductor. 

The combination of a vortex in A and a Z2 kink in A^o 
leads to an anisotropic hedgehog in n occurring where the 
vortex tube meets the interface with the topological in- 
sulator. We can think of this defect as a potential well 



VI. DEGENERACY SPLITTING 

With the Majorana wavefunctions in hand, we can use 
these expressions to calculate the splitting of the ground 
state degeneracy as a function of Majorana separation. 
This energy splitting has been calculated by a variety 
of means in 2D Majorana systems. We will generalize 
to 3D the method of Ref. |32 who adapted to 2D the 
ID Lifshitz problem^. Our calculation is very similar 
to what has already been presented in Refs. [32 and 113. 
We refer the reader to these papers for technical details. 
In what follows, we outline the main idea behind the 
calculation. 

Consider two Majorana states, and ^6, which are 
brought together from infinity. As they approach each 
other, the degenerate eigenvalue of the two fusion chan- 
nels is split by an amount ii'spiit- The new eigenfunc- 
tions of this two-Majorana state are V^j- = ^'o±i4't with 
corresponding eigenvalues E±. Particle-hole symmetry 
dictates E'sput = -E+ — = 2E+. We calculate the 
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eigenvalue of one of these states as follows^S: 



(3) 



where E is an integration region corresponding to the 
half-infinite volume in three dimensions. The integrand 
can be written in terms of total derivatives, and the lo- 
calized nature of the wavefunctions allows us to reduce 
E to the infinite 2d plane that bisects the line joining the 
two Majorana states in question. We consider two cases: 
vertical Majorana-Majorana coupling in the z-direction 
and lateral Majorana-Majorana coupling in the xy-plane. 

For the vertical coupling, we place at R-a = 
(0,0, -L/2) and *6 at R,, = (0,0,-fL/2); L is the su- 
perconductor thickness. Importantly, one of these is a 
"kink hedgehog" while the other is an "anti-kink hedge- 
hog." Equations ^ and (O lead to: 



spli 

where the prefactor is 
4V2A,, 



(^) 



10" 



sech (i/Cso) 
37 eV 



E, 



(^) 



K( ^2 / A2X g(-4/A") 



[if(-e|,/A2)-i;(-eyA2)] 



(4) 
(5) 



(6) 



While E^Q^-* is seemingly a large energy scale, the sech^ 
factor stemming from the in-plane localization makes the 
degeneracy splitting due to coupling in the z-direction 
completely negligible compared to what we calculate 
next. 

For the lateral coupling within the same interface, we 
consider two Majoranas at Ra = (— i?/2, 0, — E/2) 
and at Rh = (+i?/2,0, -L/2); R is the in-plane dis- 
tance between hedgehogs. Unlike the case of vertical cou- 
pling, the two wavefunctions here are both kinks. Equa- 
tions (ID) and ^ lead to: 



e'^^\l) con [ep R / hv p + a) _jt_ 



(7) 



where the L-dependent prefactor is 



E, 



^ 86ptanh(L/e3o)[l + (6j./A)^]-i/4 
^ ^ - [i?(-eyA2)-i^(-eyA2)] 

£;^''^^(1000A) « 10.3 meV (9) 



and the phase shift is a = (1/2) arctanei?/A. At the 
fields of interest to us, E^^f^{R, L) is much larger than 
^spiit(^)' '^^ neglect the latter. The attenuation 

and period of oscillation of E^^f^^{R, L) depends on the 
separation between Majoranas, R, which is in turn de- 
termined by the lattice spacing of the Abrikosov vor- 
tex lattice. For a triangular lattice, this spacing is re- 
lated to the field by i? ~ ^ b-!^% ' envelope 



of the energy splitting varies with field as E^^f^^^{B) cx 



2h/2e 



1.25 Tesla. This 



V'B/Boe-V-Bo/S^ where Bq 

sets the lower bound of the temperature window and 
is clearly field-dependent: Tl[B) ~ E^^i^^{B)/kB- See 
Fig. [3^. We next determine the upper bound of the tem- 
perature window. 



VII. ENTROPY 

To compute the cooling effect we need to understand 
all appreciable contributions to the total entropy as a 
fimction of temperature and field: S{T,B). These can 
be classified into phonon (Sph) and vortex (Sy) contri- 
butions with the latter being composed of several pieces 
(electronic contributions at the temperatures of interest 
are negligible because A, A^o ^ ksT). Thus, the total 
entropy is given by 



S{T,B)^SphiT) + SyiT,B) 
The phonon entropy is standard: 

47r4y / T 



SphiT) 



5K,„ 



e 



D 



(10) 



(11) 



where &p) ~ 182 K is the Debye temperature of the 

parent compound^ and Vuc ~ 426 A is the volume 
of the unit cell^^. Within the superconductor material 
Cu3;Bi2Se3, recent specific heat data3^ has yielded a De- 
bye temperature oi Qd ^ 120 K. 

The vortex entropy will have several contributions. 
The most important piece of Sy is due to the non- 
Abelian anyons (5„a); this is what drives the dramatic 
low-temperature cooling effect. For a large number of 
vortices (Ny) this is simply^ 



Sna{B) ^ Ny{B)2kB \ogV2. 



(12) 



The non-Abelian anyon entropy only depends on B, not 
T, and this is through its dependence on the number of 
vortices: Sna{B) cx Ny{B) with Ny{B) k, BA/4)q where 
(fn^ ~ h/2e\s the flux quantum and A is the sample area 
perpendicular to the magnetic fleld. This linear approx- 
imation is justified for B 3> B^i- 

The contribution to the vortex entropy from more con- 
ventional sources will depend on the vortex density, and 
thus the magnetic field. An isolated vortex line will con- 
tribute two types of entropy. First, there are subgap 
bound states that are localized to the core but extended 
along the z-direction. These are usually called CdGM 
excitations after Caroli, de Gennes, and Matricon^^. Sec- 
ond, a single vortex line will have excitations analogous 
to those of a fluctuating string^!. For our materials and 
parameter regimes, this second type of fluctuation turns 
out to contribute negligibly to the total entropy. 

In addition to this single-vortex physics, collective ef- 
fects can manifest when the magnetic field is increased 
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even slightly above Bd leading to the formation of an 
Abrikosov vortex lattice. A dense array of vortices can 
have collective modes of the same two types as described 
above. First, there are collective CdGM excitations^, 
and second there are modes corresponding to fluctu- 
ating elastic mediai^. Note, however, that the term 
"dense" must be understood with respect to the appro- 
priate length scale. Define R as the lateral vortex-vortex 
distance, A as the penetration depth, and ^ as the su- 
perconducting coherence length. Collective modes of the 
vortex lattice appear when the magnetic field is such that 
R < X. In contrast, collective CdGM excitations only ap- 
pear at the much higher density R ^ ^. For our device, 
Bel ~ 1.3 mT and Bc2 = 1-7 T, giving a rather large 
range R X (or, equivalently, Bci <C -B ^ Bc2) in 
the dilute limit with respect to CdGM excitations, but 
the dense collective-mode limit of vortex fluctuations. 

In such a regime, the CdGM entropy is given by Ny 
times the "isolated" vortex line contribution. 



(rr. m Ny{B)kB2kFL r~9~-g/kuT /^J^ 

5cdGM(r,i?) = -= J^e (13) 



5/2 

(14) 



while the vortex lattice entropy takes the forntiii 



Syi{T,B) = fcB^C(5/2) 



3/2 



where g = 2^:^ is the CdGM mini-gap, k = h/2m « 2 x 

10^ ''A'^/s is the quantum of circulation, ^ is the Riemann- 
Zeta function, and V = LA is the superconductor sample 
volume. This form of Syi, which is proportional to _B^^/^, 
is only valid in the parameter regimes of interest to us. 
Eventually, at very low fields on the order of Bd, S^i 
must of course vanish as B decreases^. 
Thus, the total vortex entropy is given by 

SyiT, B) = S,,a{B) + 5cdGM(T, B) + S,a{T, B) (15) 

When added to the phonon entropy, this yields an ap- 
proximate analytic expression for the total system en- 
tropy as a function of T and B. Using these expressions 
we calculate the central quantity ^ = — [ds/drjfl 
scribed in the next section. 



VIII. COOLING 

To understand the cooling effect, consider the small 
change in entropy for a system depending on temperature 
and field: 



dS 



dS_ 

dT 



dT 



dS_ 
dB 



dB 



(16) 



For an isentropic process {dS = 0), the system's tempera- 
ture changes in response to a small field change according 
to: 



dT 
dB 



{dS/dB)j 
{dS/dT)i 



Without 
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T(K) 



(17) 



FIG. 2: (color online) Quantitative illustration of cooling via 
3D non-Abelian anyons. The superconductor layer thickness 
is set to L = 1000 A. The orange and purple dashed lines are 
constant field cuts {B = 25-Bci ~ 33 mT) of ^ correspond- 
ing to the orange and purple dashed lines in Fig. [3] with (a) 
and without (b) anyons. Notice that Th <i Tu, where Th is 
defined by the temperature below which |^ > in a system 
without 3D Majoranas, and Tu, the upper bound of the cool- 
ing window, is defined as the temperature below which the 
non-Abelian entropy strongly enhances the cooling effect in a 
system with 3D Majoranas. 



When this quantity is negative it represents a decrease 
of system's temperature as the field is adiabatically in- 
creased: cooling. In contrast, a positive sign indicates 
heating. Importantly, the sign and magnitude of ^ de- 
pend on T and B. Since ^ is always positive, to find 
cooling we require a parameter regime in which ^ is pos- 
itive. This will occur when Sna{B) dominates the total 
entropy. 

In Fig. [5] we show the temperature dependence, both 
with and without 3D anyons in the system, for a fixed 
value of field: B = 25Bci ^ Bd- At high tempera- 
tures dT /dB is very similar in both cases, but below a 
certain temperature, Th, the system without 3D anyons 
will experience heating while the system with 3D anyons 
will experience cooling. This qualitative difference is the 
harbinger of 3D non-Abelian anyons. 

To reiterate: weak cooling may occur at high tempera- 
tures with or without 3D anyons, but in a specific region 
of the iJT-plane, depicted in dark blue in Fig[3^, a system 
with 3D anyons will exhibit strong cooling while a system 
without 3D anyons will experience heating (as shown in 
Fig Wp). For a system with 3D anyons, the transition 
between the strong cooling and weak cooling regimes is 
defined by the extremum of the curve dT /dB which is 
most apparent by examining the purple dashed curve in 
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FIG. 3: (color online) Comparison of the temperature evolu- 
tion with (a) and without (b) 3D non-Abelian anyons in the 
field-temperature plane. The superconductor layer thickness 
is set to L = fOOO A. The sign and magnitude of |§(B,T) 
determines the nature of the cooling or heating effect as the 
magnetic field is adiabatically increased. When 3D Majoranas 
are present (a), a sizable temperature window Tl < T < Tu 
at intermediate fields (darker blue region labelled "Strong 
Cooling") will exhibit a strong cooling effect. Without 3D 
Majoranas (b), the same parameter regime will demonstrate 
a heating effect. The yellow region below Tl{B) is where the 
ground state degeneracy is split by Majorana-Majorana tun- 
neling; oscillations are not shown for clarity. This interesting 
sector may include collective anyon excitations generalizing 
some ideas in 2D^, but these are not pertinent to this paper. 

Fig.m 

To understand the origin of the minor heating with- 
out anyons note that, in the regimes relevant to our pro- 
posal, all contributions to the entropy are non-decreasing 
functions of B except S^i. This part of the entropy de- 
creases with field and thus dS/dB could become negative 
if Svi were large enough compared to all the other entropy 
sources. Since dS/dT is always positive, the possibility 
of negative dS/dB means dT/dB = —^§7^ could be 
driven positive leading to the observation of heating with 
an adiabatic increase of the field. Indeed, in the absence 
of 3D non-abelian anyons this is exactly what happens 
at low temperatures as depicted in Fig.[3)D. However, the 
presence of 3D non-abelian anyons contributes Sna to the 
total entropy which increases sufficiently rapidly with B 
that it overwhelms the low temperature heating due to 
Syi and reverses the trend to produce a dramatic strong 
cooling effect. 

IX. CONCLUSION 

In summary, we have quantitatively estimated a 
magnetic field and temperature regime in which an 
experiment might detect 3D non-Abelian anyons in 



a heterostructure device composed of Bi2Se3 and 
Cug i2Bi2Se3. Within this regime, a system with 3D non- 
Abelian anyons will experience a decrease in temperature 
as the magnetic field is adiabatically increased. In con- 
trast, a system without 3D non-Abelian anyons under 
identical conditions will exhibit a temperature increase 
as the magnetic field is adiabatically increased. In the 
strong cooling region, the effect is rather large. For ex- 
ample, at B 25Bci and T = 0.1 K we have |g w -1.7 
K/T, which should be detectable with present technol- 
ogy 

We close by enumerating a few caveats that are spe- 
cific to Cuq. 12612863, but not to the general idea of our 
proposal. First, the numbers presented in this paper are 
based on a value of the superconducting gap inferred from 
the i3c2-derived coherence length^^ which disagrees with 
the Tc derived gap2^. Since the degeneracy splitting (i.e. 
the Tl{B) line) is exponentially sensitive to this param- 
eter, the choice is very important. We hope a direct ex- 
perimental measurement of the superconducting gap will 
be made in the future to pin down the true value of this 
important parameter. Second, we have assumed the sign 
of Aso to take opposite values in insulating Bi2Se3 versus 
Cuq. 12612863. Without the sign change, there will be no 
Majorana mode. Again, this needs to be experimentally 
checked for Cuq. 12612863. Note, however, that if later 
experimental investigations determine that either A in 
Cuq. 12612863 is much smaller than our assumption, or 
that Ago does not change sign in Cuq. 12612 803, the gen- 
eral idea of this proposal will not be invalidated but only 
its applicability to this particular superconductor; all 
qualitative conclusions will remain true for any insulator- 
superconductor system that satisfies the following condi- 
tions: (i) the superconductor must have a relatively large 
s-wave gap and be strongly type II; (ii) the band gap 
must take opposite signs in the insulating and supercon- 
ducting regions. These are relatively simple conditions 
to fulfill. We chose to examine a 6i2Se3-based structure 
because this material has emerged as an archtypc topo- 
logical insulator which is being independently studied by 
many different research institutions. Furthermore, the 
possibility of creating topological insulator regions and 
superconducting regions simply by p-doping or n-doping 
the same parent compound makes 612803 a very attrac- 
tive system for building heterostructure devices in the 
future. 
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